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ABSTRACT 

Theoretical modeling of the driving processes of solar-like oscillations is a powerful way of understanding the properties of the 
convective zones of solar-type stars. In this framework, the description of the temporal correlation between turbulent eddies is an 
essential ingredient to model mode amplitudes. However, there is a debate between a Gaussian or Lorentzian description of the eddy- 
time correlation function (Samadi et al. 2003, Chaplin et al. 2005). Indeed, a Gaussian description reproduces the low-frequency 
shape of the mode amplitude for the Sun, but is unsatisfactory from a theoretical point of view (Houdek, 2009) and leads to other 
disagreements with observations (Samadi et al., 2007). These are solved by using a Lorentzian description, but there the low-frequency 
shape of the solar observations is not correctly reproduced. We reconcile the two descriptions by adopting the sweeping approximation, 
which consists in assuming that the eddy-time-correlation function is dominated by the advection of eddies, in the inertial range, by 
energy-bearing eddies. Using a Lorentzian function together with a cut-off frequency derived from the sweeping assumption allows 
us to reproduce the low-frequency shape of the observations. This result also constitutes a validation of the sweeping assumption for 
highly turbulent flows as in the solar case. 
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1. Introduction 

Excitation of solar-like oscillations is attributed to turbulent mo- 
tions that excite p modes (for a recent review, see Samadi 2009]). 
Their amplitudes result from a balance between excitation and 
damping and crucially depend on the way the eddi es are tempo- 
rally correlated as show i i for so l ar p and ^ modes ( Samadi et al.l 



drodynamical simulations of the outer layers of the Sun show 
that at the length-scales close to those of the energ y-bearing ed 



dies (around 1 Mm), yj; is a L orentzian function dSamadi et al 



l2003bt IBelkacem etaP '2009b': 'Appourchaux et al.l 1201 Ol) . for 
main-sequence sta rs (|Samadi et al. 20 10b, a ), for red giants 
(iDupret et al.l2009h . or for massive stars (IBelkacem et al.l2009ai 
120101) . Hence, the improvement of our understanding and mod- 
eling of the temporal correlation of turbulent eddies, hereafter 
denoted in the Fourier domain as Xk{<^), is fundamental to infer 
turbulent properties in stellar convection zones. 

There are two ways to compute the eddy-time correlation 
function. A direct computation from 3D numerical simula- 
tions is possi ble and was pe rformed by Samadi et al. (2008^. 
Nevertheless. ISamadil (l2009l) pointed out that the results depend 
on the spatial resolution, and therefore dedicated high-resolution 
3D numerical simulations are required. This then becomes an 
important limitation when computing mode amplitudes for a 
large number of stars, preventing us from applying statistical 
astereosismology. 

The second way to compute Xk consists in using appro- 
priate analytical descriptions. Most of the theoretical formu- 
lations of mode excitation explicit ly or implicitly assume a 
Gaussian functional form for Xkiio) (Goldreich & Keelev 1977t 

till 



Dolginov & MusUmovll 19841 iGoldreich et a 



199ilSamadi et al.l2001l;IChapUn et al.ll2005l) . However, 3D hy 
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l2003at IBelkacem et al.ll2009bl) . As pointed out bv lChapUn et al 
(l2005h . a Lorentzian is also a result predicted for the largest, 
most-energetic eddies by the time-dependent mixing-length for- 
mulation derived by|Gpugh (197j). Therefore, th ere is numer- 
ical, t heoretical, and also observational evidence dSamadi et al.l 
120071) that;^f^: is Lorentzian. 

Ho wever, IChaplin et al. l (l2005l) . [Samadil(l2009l) . and lHoudekl 
(120091) found that a Lorentzian Xk, when used with a mixing- 
length description of the whole convection zone, results in a se- 
vere over-estimate for the low-frequency modes. In this case, 
low-frequency modes (v < 2mHz) are excited deep in the so- 
lar convective region by large-scale eddies that g ive a substan- 
tial fra ctio n of the energy injected to the modes. IChaplin et al.l 
(120051) and lSamadil (l2009h then suggested that most contributing 
eddies situated deep in the Sun have a Xk more Gaussian than 
Lorentzian because at a fixed frequency, a Gsaissianxk decreases 
more rapidly with depth. 

We therefore propose a refined description of the eddy-time 
coiTelation function based on the sweeping approximation to 
overcome this issue. This consists in assuming that the temporal 
correlation of the eddies, in the inertial subrange, is dominated 
by the advection by energy-bearin g eddies. This assumption was 
first prop osed by TennekesI d 19751). and was subse quently inves- 
tigated bv lKanedal(ll993h and lKaneda et al] (Il999l) . In this letter, 
we demonstrate that the low-frequency shape of the observed en- 
ergy injection rates into the solar modes is very sensitive to this 
assumption and more precisely to the Eulerian microscale, de- 
fined as the curvature of the time-correlation function at the ori- 
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gin. Hence, modeling of the solar p-mode amplitudes is shown 
to constitute an efficient test for temporal properties in highly 
turbulent flows. 

The paper is organized as follows: Section |2] defines the 
eddy-time correlation function. In Sect. [3] we propose a short- 
time expansion of the eddy-time correlation function. The use of 
the Eulerian microscale as a cut-off is introduced in the compu- 
tation of solar p mode amplitudes and the result is compared to 
the observations in Sect. |4] Finally, Sect. |5] is dedicated to con- 
clusions and discussions. 



2. The Eulerian eddy time-correlation function 

For a turbulent fluid, one defines the Eulerian eddy time- 
correlation function as 



{u{x + r,t + T)- u(x, 0) = I S(k, t, r) e'*'^ d^Jt , 



(1) 



where u is the Eulerian turbulent velocity field, x and t the space 
and time position of the fluid element, k the wave number vec- 
tor, T the time-correlation length, and r the space-correlation 
length. The function & in the RHS of Eq. ([T]i represents the time- 
correlation function associated with an eddy of wave-number k. 

We assume an isotropic and stationary turbulence, accord- 
ingly £ is only a function of k and t. The quantity &{k, r) is 
related to the turbulent energy spectrum according to 



S(k,T) 



E(k, t) 
Ink^ 



(2) 



where E{k,T) is the turbulent kinetic energy spectrum whose 
temporal Fourier transform is 



E{k 



Eik, t) e"^^ dr , 



(3) 



where cj is the eddy frequency, an d E(k,a)) is w ritten following 
an approximated form proposed bv lSteinI (Il967h 



E{k, oS) - E(k)xk{^) with 



f 



Xk{c^)du^\, (4) 



where Xk((Jj) is the frequency component of E{k, oS). In other 
words, represents - in the frequency domain - the tempo- 
ral correlation between eddies of wave-number k. 

As already discussed in Sect.[Tl theoretical and observational 
evidence show that;^f^(w) is Lorentzian, i.e. 



1 



1 



^ii^k 1 -I- (_u)la)k) 



(5) 



where lOk is by defi nition the width at half m aximum of Xk{(J^)- 
In the framework of ISamadi & Goupill (1200 ih 's formalism, this 
latter quantity is evaluated as: 



u)k - kuk with 



f-2k 



E{k)dk, 



(6) 



where E{k) is defined by Eq. However, in the high-frequency 
regime (a> » ojk), corresponding to the short-time correlation 
(t X 0), the situation is less clear We next investigate short-time 
correlations (r x 0). 



3. The sweeping assumption for the Eulerian 
time-correlation function 

3. 1 . Short-time expansion of tine eddy-time correlation 
function 

The function &{k, t, r) (Eq. ([TJ) can be expanded for short-time 
scales, in the inertial sub-range {i.e for k > ko and k < k^, where 
ko is the wave number of energy-bearing eddies and k^ is the 
wave-number of viscous dissipation), using the Navier-Stokes 
equations and the sweeping assumption, as (see lKanedairi993L 
for a derivation) 



6{k, r) = 6{k, T = 0) |l - ak\T\ - 1 (ojet)- + . . . j 



(7) 



where the characteristic frequency ak is defined by the relation 



= --—(«■«)= rQ'*fi(A:,T = 0)d^jt 
2 df J 



(8) 



with 6 the dissipation rate of energy. Hence, ak is the typical 
frequency of energy dissipation at the scale k. It can be esti- 
mated by assuming that a large fraction of the kinetic energy 
of ed dies is lost within one turnover time ( iTennekes & Lumlevi 
Hence, ak is approximated by the turn-over frequency 
Uk - kuk - LOk (see Eq. (|6]l). 

The second characteristic frequency, a)E{k), is the curvature 
of the correlation function near the origin (iKanedail993i) . and is 
defined by 



ue -kuQ. 



(9) 



The associated characteristic time TE{k) - Inu^^ is also 

referr ed to as the Eulerian micro-scal43 (ITennekes & Lumlevi 
1 19721) . An approximate expression for (x)E{k) can be obtained by 
assuming the random sweeping effect of large eddies on small 
eddies. This assumption consists in assuming that the velocity 
field u{k) associated with an eddy of wave-number k lying in 
the inertial-subrange {i.e. large k compared to ^o) is advected by 
the energy-bearing eddies with velocity hq {i.e. of wave-number 
^:o)- This time-scale is obtained by assuming uniform density, 
which is valid in the Sun for k > ko {i.e. in the inertial sub- 
range) since the density scale-height approximately equals the 
length-scale of energy-bearing eddies ((dlnp/dr)"' Inlko). 
It also assumes the quasi-normal approximation ( see iKanedal 
ll993l:lKaneda et al.lll999HRubinstein & Zhoul2002l for details). 
The Eulerian micro-scale then corresponds to the timescale over 
which the energy-bearing eddies of velocity uq advect eddies 
of size Ink'^ . It is the lowest time-scale (highest frequency) on 
which those eddies can be advected. 



3.2. Eulerian time micro-scale as a cut-off frequency 

The issue is now to estimate to what extent cl>e can be considered 
as a cut-off frequency, i.e. that there is a sharp change in the slope 
of Xk at high co. 

To this end, we first remark that the zero-th- and first-order 
terms in Eq. (|7]i are consistent with an exponential decrease of 
width Ok {i.e. a Lorentzian in the frequency domain of width 
wa, Eq. (|5]l) for small r. In contrast, the zero-th-order term to- 
gether with the second order term in Eq. (Q are consistent with 



' It is tlie time equivalent of the Taylor micro-scale, which corre- 
sponds to the largest scale at which viscosity affects the dynamic of 
eddies. 
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Fig. 1. Schematic time-correlation (xk) versus normalized eddy 
frequency (cj/oje) at ^ = 5^o ('-e. in the inertial subrange such 
as kd :» k > ko ), where ko - 6.28 x lO^^'m"'. Note that the 
value of ^0 does not influence the result. The solid line (resp. 
dashed triple dot line) correponds to the Lorentizan functional 
form of Xk for o) < a>E (resp. ai < coe). The dashed line corre- 



sponds to a Gaussian modeling (xk °^ e 



-(oj/ojeY 



) of characteristic 



frequency oie- (we numerically verified that e"('"/'"£)" is a good 
approximation of Eq. (|7]i, see also Sect. I3.2l i. We stress that the 
sharp decrease the functional form given by Eq. (|7]), in the tem- 
poral Fourier domain, then justifies to consider oje as a cut-off 
frequency. In other words, Xk is computed according to Eq. (fT2l) . 



a Gaussian behavior of width te- Hence, the relative importance 
of those two regimes depends on the relative magnitude of the 
second and third terms in Eq. (|7|. Let us define the ratio (??) of 
the first to the second order term in the expansion of £ (Eq. ^) 



9? = 2 (w£t)- 



CL>k 

a>E 



(10) 



To evaluate this ratio, we compare the typical frequencies ojk and 
oje using Eq. ^ together with Eq. Adopting a Kolmogorov 
spectrum E{k) = Ck e^^^ k^^^^, with Ck the Kolmogorov univer- 
sal constant (close to 1.72), we have u\— fiu^ {k/koy^^^ , where 
/3 = 0.555. Hence 



coe 
a>k 



uq 

Uk 



.-1/2 



1/3 



(11) 



From Eq. ( fTTT i we conclude that for k :i> ko (i.e. at small scale) 
we have w^/w^ x /3'^^~ (k/koY^^ » 1, then te «: And for 
k X ko (i.e. at large scale), we have w^/wa. = yS"'^^ ~ 1.4. Hence, 
we always are in the situation where oie > oJk- 

From Eq. (fTOl i. it immediately follows that for u) > oje the 
second order term dominates over the first order one in Eq. (|7]i, 
at all length-scales. We then conclude that for frequencies near 
the micro-scale frequency {co > cje), the eddy-time correlation 
function behaves as a Gaussian function (e"("/"t)') instead of a 
Lorentzian function, resulting in a sharp decrease with to (see 
Fig. [TJ. Hence, the contributions for w > are negligible and 
the temporal correlation is computed as follows 



1 



1 + {to/aJkY 



if w < (aJe 



(12) 



if oj > oje ■ 



4. Computation of the p-mode energy injection 
rates 

4.1. Computation of the energy injection rate 

The formalism we used to compute ex citation r ates o f ra- 
dial modes was de velop ed by .S am adi & Goupill (1200 Ih and 
ISamadi et al.l (l2005h (see I Samadi 2009L for a thorough discus- 
sion) 

For a radial mode of frequency wo = 27r vq, the excitation 
rate (or equivalently, the energy injection rate), P, mostly arises 
from the Reynolds str esses and can be written as (see Eq. (21) of 
iBelkacem etal]|2008h 



Sk 



21 Jo 

EHk) 
k^ 



Po 



dr 



Sk dk 



X+o. 
oo 



dm (13) 



(14) 



where is the radial component of the fluid displacement eigen- 
function (^), m is the local mass, po the mean density, the 
mode angular frequency, / the mode inertia, Sk the source func- 
tion, E{k) the spatial kinetic energy spectrum, Xk the eddy-time 
correlation function, and k the wave-number 

The rate (P) at which energy is injected into a mode is com- 
puted according to Eq. ( fT3] l. In this letter, we consider two theo- 
retical models, namely: 



An analytical approach: the ID calibrated 
ture model used for these computations 
with the stellar evolution code CESAM2k 



solar struc- 
is ob tained 
(iMorell 



iMorel & Lebretonll2008l) . The atmosphere is computed as- 
suming an Eddington grey atmosphere. Convection is in- 
cluded accordi ng to a Bohm-Viten se mixing-length (MLT) 
formalism (see lSamadi et al.ll2006l for details), from which 
the convective velocity is computed. The mixing-length pa- 
rameter a is adjusted in a way that the model reproduces 
the solar radius and the solar luminosity at the solar age. 
This calibration gives a = 1.90, with an helium mass 
fraction of 0.245, and a chemical composition following 
iGrevesse & Noelsl (Il993h . The equilibrium model also in- 
cludes turbulent pressure. 
- A semi-analytical approach: calculation of the mode ex- 
citation rates is per formed essentially in the manner of 
ISamadi et al.l ( l2008allbl) . All required quantities, except the 
eddy-time correlation function, the mode eigenfunctions (^,.) 
and mode inertia (/), are directly obtained from a 3D simu- 
lation of the outer layers of the Sun (see [Samadi., 20091 for 
details on the numerical simulation). 

In both cases, the eigenfrequencies and eigenfunctions 
are computed with the ad iabatic pulsation code ADIPLS 
dChristensen-Dals gaardl2008h . We stress again that in both cases 
Xk is implemented as an analytical function. 

4.2. Results on mode amplitudes 

When the frequency range of Xk is extended toward infinity, 
computation of P according to Eq. (fT3T l and Eq. (fT4l i fails to re- 
produce the observations, in particular the low-frequency shape. 
In order to illustrate this issue, we have computed the solar 
model excitation rates, using the solar global model described 
in Sect. 14.11 The turbulent kinetic energy spectrum {E{k)) is as- 
sumed to be a Kolmogorov spe ctrum to be co nsistent with the 
derivation of te as proposed bv lKanedal (Il993h . In addition, the 
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eddy-time correlation function is supposed to be Lorentzian as 
describ ed by Eq.jSl l for all oj > 0. In agreement with the re- 
sults of lChaphneta l. (2005) andlSamadi (2009), it results in an 
over-estimate of the excitation rates at low frequency (see Fig.|2] 
top). 

In contrast, by assuming that the time-dynamic of eddies 
in the Eulerian point of view is dominated by the sweeping, 
the Eulerian time micro-scale arises as a cut-off frequency (see 
Sect. l3.2l i. Hence, yi-(oj) is modeled following Eq. (|5]) for w < oje 
andxk(<^) = elsewhere. Using this procedure to model Xk{<^) 
{i.e. by introducing coe as a cut-off frequency) permits us to re- 
produce the low-frequency (v < 3mHz) shape of the mode exci- 
tation rates as observed by the GONG network (see Fig.|2]i. This 
is explained as follows: for large-scale eddies near A:^', situated 
deep in the convective region, the cut-off frequency oje is close 
to ail; as shown by Eq. (fTTI) . As a consequence, the frequency 
range over which ji't is integrated in Eq. ( fT4l) is limited, resulting 
in lower injection rates into the modes. 

Note that the absolute values of mode excitation rates are not 
reproduced by using a mixing-length description of convection, 
this is in agreement with Samadi (2009), and arises because that 
it underestimates the convective velocities as well as convective 
length-scales. It then explains the differences between the com- 
putation of mode excitation rates using the MET and the 3D nu- 
merical simulations (Fig.|2l). 

5. Conclusion and discussion 

By using a short-time analysis and the sweeping assumption, we 
have shown that there is a frequency {cje the micro-scale fre- 
quency) beyond the temporal correlation Xk sharply decreases 
with frequency. Including this frequency as a cut-off in our mod- 
eling of Xk and assuming a Lorentzian shape we are able to 
reproduce the observed low-frequency (y < 3 mHz) excitation 
rates. 

These results then re-conciUate the theoretical and observa- 
tional evidence that the frequency dependence of the eddy-time 
correlation may be Lorentzian in the whole solar convective re- 
gion down to the cut-off frequency coe- Finally, it also repre- 
sents a validation of the sweeping assumption in highly turbulent 
flows. 

We note, however, that one must remove several theoretical 
shortcomings to go further. For instance, a rigourous treatment 
of the energy-bearing eddies is needed. The short-time analysis 
and the computation of the Eulerian miscro-scale must be recon- 
sidered by including the effect of buoyancy that mainly affects 
large scales. Furthermore, some discrepancies remain at high- 
frequency (v > 3 mHz), and to go beyond these one ha s to re- 
move the separation of scales assumption (see Belkace m et alJ 
120081 for a dedicated discussion) and include non-adiabatic ef- 
fects. 

Eventually, we note that the modeling of amplitudes under 
the sweeping assumption is to be extended. In particular, the 
investigation of the effect of the sweeping assumption on solar 
gravity mode amplitudes is desirable. 
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